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THE DISTRIBUTION OF PERMUTATION MATRIX ENTRIES 
UNDER RANDOMIZED BASIS 

BENJAMIN TSOU 


Abstract. We study the distribution of entries of a random permutation 
matrix under a “randomized basis,” i.e. we conjugate the random permuta¬ 
tion matrix by an independent random orthogonal matrix drawn from Haar 
measure. It is shown that under certain conditions, the linear combination of 
entries of a random permutation matrix under a “randomized basis” converges 
to a sum of independent variables sY + Z where Y is Poisson distributed, Z 
is normally distributed, and s is a constant. 


1. Introduction 

Traditionally, random matrix theory has focused primarily on the eigenvalue dis¬ 
tributions of matrices drawn from various measures (see e.g. Q for a quick survey or 
[l| for an introduction to the held). However, there has also been interest in study¬ 
ing statistics related to entries of random matrices. For example, the study of the 
entries of a Haar-distributed orthogonal matrix is a classical subject originated by 
Borel 0 who showed that the normalized hrst coordinate of a random vector chosen 
from the n-dimensional sphere converges in distribution to a standard normal. Di- 
aconis and Freedman [l^ strengthened this to obtain total variation convergence of 
the hrst k coordinates to independent standard normals when k = o{n). Jiang [l^ 
then further strengthened this result considerably to show that if = o{^/n) and 
qn = o{y/n), then the total variation distance between the joint distribution of the 
entries of the upper left x Qn block of a Haar-distributed orthogonal nxn matrix 
and PnQn independent standard normals converges to 0. Moreover, he proved this 
result is sharp. 

Matrix entries of random permutation matrices have also been studied. Hoeffd- 
ing’s combinatorial central limit theorem [l5j gives a result on linear combinations 
of random permutation matrix entries. Let A„ be a sequence of n x n real matrices 
such that Tr(A„A^) = n. Hoeffding shows that under certain conditions on A„, 
if P is a random nxn permutation matrix, then Tr(A„P) converges weakly to a 
standard normal random variable as n —^ oo. 

In this article, we study the distribution of entries of a random permutation 
matrix in a basis free way. Individual entries of a permutation matrix are either 
0 or 1, but if we conjugate by a random orthogonal matrix, we can look at the 
distribution of entries under a “randomized basis.” To formalize this, let P be an nx 
n permutation matrix drawn uniformly from the symmetric group ©„ and let M be 
a Haar-distributed nxn orthogonal matrix independent of P. (With little confusion, 
we will often conflate a permutation matrix and its corresponding permutation in 
©„). Let An be a sequence of n x n real matrices such that Tr{AnA^) = n. Then 
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it is shown in Theorem 13.21 below that under certain conditions on A„, 

Tr{AnMPM^) A sY + Z 

as n —>■ oo where V and Z are independent random variables such that Y is Poisson 
distributed, Z is normally distributed, and s is a constant. Whereas Hoeffding’s 
combinatorial CLT shows convergence for linear combinations of random permuta¬ 
tion matrix entries, Theorem l3.2l can be interpreted as a distributional convergence 
for linear combinations of random permutation matrix entries under a “randomized 
basis.” 

We briefly review some other related results. Let A be an n x n real (nonran¬ 
dom) matrix such that Tr{AA'^) = n. D’Aristotile, Diaconis, and Newman 
have shown using characteristic function methods that if M is an n x n random 
orthogonal matrix, then Tr(AM) converges in distribution to a standard normal 
random variable as n —>■ oo uniformly in A. Recently, Meckes 0 has improved this 
result. She shows that the total variation distance between Tr{AM) and a stan¬ 
dard normal random variable is bounded by and that this rate is sharp up to a 
constant. In Q, Chatterjee and Meckes obtain bounds on the Wasserstein distance 
between the multivariate distributions {Tr{AiM),Tr{A 2 M), ...,Tr{AkM)) and a 
Gaussian random vector. Their techniques involve generalizing Stein’s method of 
exchangeable pairs. 

2. Limiting distribution of a single entry 

In this section, we start by showing that a single scaled matrix entry of MPM'^ 
converges in distribution to a standard normal distribution. Later sections will 
consider joint distributions of the matrix entries. 

Note that the distribution of where P is random can be thought of as a 

mixture of distributions of MPM"^ where P is a hxed permutation matrix. Thus, 
we first assume P is fixed. If P is similar to another permutation matrix Q via an 
orthogonal transformation V, i.e. P = VQV^, then MPM"’' = {MV)Q{MVY’ = 
MQM'^. Thus, the distribution of MPM"'^ only depends on the cycle type of the 
permutation matrix P. 

First, we consider the case where the permutation consists of the single n-cycle 
(l,2,...,n). Let Cn be the corresponding permutation matrix. Then 

n 

{MCnM^)ab = E 

(where abusing notation slightly, := i) By symmetry, it is clear that 

] = 1/n. Moreover, \/n{Ma,i, (Ai, A 2 ) where Xi and X 2 are i.i.d 

A/’(0,1). (The convergence is actually in total variation, see e.g. 0,0). Therefore, 

n 

should have variance on the order of I/n and it is natural to consider 

i—1 

the scaled random variable \/n{MCnM^)ab- 

To obtain the limiting distribution of ^/n{MCnM'^)ab (as well as the limiting 
joint distribution of multiple entries derived later on), the key idea will be to use a 
martingale central limit theorem for dependent random variables proved by McLeish 
0 . Specifically, we use his Corollary (2.13). (See 0 for a book-length treatment 
on martingale limit theorems.) Let {Xn,i] 1 < * < n} be a triangular array of 
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n n 

random variables. Define ^ Let 

2=1 2—1 

= o’(^n,i, ^n, 2 , We say that Xn^i is a martingale difference array if 

V.{Xn^i\J-n,i-i) = 0. Then we have the following: 

Theorem 2.1 (McLeish). Suppose X„^i is a martingale difference array normalized 

n „ 

by its variance satisfying the Lindeberg condition lim N / X^ ^dP = 0 

for all e > 0 and the condition lim sup < 1- Then Sr. 4 Ar(0,1). 

n—>-oo . . . 


The following lemma will also be useful in the sequel: 


Lemma 2.1. Let the random variable Xn = where ki + ... + kr = 2k 

and Mi are distinct entries of the nx n Haar distributed orthogonal matrix. Define 

Y = where Zi are i.i.d. standard normals. Then Xn —^ Y. 


Proof. By [l^ , the joint distribution of any k fixed entries (scaled by y/n) of a Haar 
distributed orthogonal matrix converges in total variation to the joint distribution 
of k i.i.d. standard normals. In particular, by the continuous mapping theorem, we 
see that Xn 4 Y. To get convergence, we need uniform integrability, which is 
implied by limsupE[X^] < oo. By the Cauchy-Schwarz inequality, it is sufficient to 

n—foo 

show that limsupE[n^^M4] < oo. But Mf has a Beta distribution with parameters 

n—¥oc) 

(1/2, (n — l)/2) (see e.g. [1^ pp. 145]) and hence E[n^^M4] = 0(1). 

□ 


Remark 2.1. Lemma [Q gives rather crude asymptotics for the expectation of 
products of Haar distributed matrix entries. Collins and Sniady @ have devel¬ 
oped machinery using Weingarten functions to compute the integrals of polynomial 
functions of Haar distributed orthogonal, unitary, and symplectic entries for any 
dimension n. (Please refer to @ for the definition of Weingarten functions. We will 
only require the asymptotics of Weingarten functions mentioned below.) 

For the orthogonal group, Collins and Sniady’s formula says that the integral 

over Haar measure / Mirjr...Mi.^^j.^^dM is the sum of Weingarten functions 

JM& 0 {n) 

W gfm., n) over pair partitions m and n of the set {1,..., 2r} such that irn(2fc-i) = *m(2fc) 
and j„(2fe-i) =jn(2fc) where m = {{m(l), m(2)},..., {m(2r - 1), m(2r)}} and n = 
{{n(l), n(2)},..., {n(2r —1), n(2r)}}. In other words, m represents pairing up matrix 
entries in the same row and n represents pairing up entries in the same column. In 
particular, the expectation is non-zero only if there are an even number of entries 
in each row and column. 

In 0 , Collins and Matsumoto obtain a simpler formula for the orthogonal Wein¬ 
garten function that significantly reduces the complexity involved in the compu¬ 
tation of Weingarten formulas and show that the integrals are always rational 
functions of the dimension n. They also provide asymptotics for the Weingarten 
function: 

Wg(m,n) = 0(l/n-2"+'(’"’'^)) 

Here, the length Z(m, n) is defined as follows: consider the graph with vertex set 
1,2,..., 2r and edge set that consists of {m(2A: — 1), m(2fc)}, {n(2fc — 1), n(2fe)} where 


4 


BENJAMIN TSOU 


1 < fc < r. Then Z(m,n) is the number of components in this graph. Thus, finding 
the asyiuptotics of the integral / is equivalent to finding 

JM^O(n) 

the (admissible) pairings m and n that maximize the number of components in 
the graph. Q also contains a useful table of values for the orthogonal Weingarten 
functions of small dimension. 


For our purposes, it will be simpler to use the stronger Lyapunov’s condition 
in place of the Lindeberg condition in McLeish’s theorem. For 5 > 0, Lyapunov’s 
condition says that 


lim 

n—¥cci 








|2+J 


= 0 . 


( 2 . 1 ) 


We apply this now to show 


Lemma 2.2. 

n 

v^E Af(0,l) 

i=r 

Proof. Let = v^ALa.i-^b.i+r- By sign syiumetry, = 0. Note that the 

random vector •■•Wn,*-i) is only a function of (Ma,i,MQ_ 2 , 

Since does not share a column with any of these ma¬ 

trix entries, its sign can be reversed independently while preserving Haar measure. 
Thus, 

{Xn^lj X.12^2^ —1; ibffoy+l ) — {Xn^lj X.it^2j —1; il^fc,i+r) 
and therefore 


= y/n¥.{Ma,iMb,i+l\(j{Xn^l,Xn^2, ■■■Xn,i-l)) = 0 . 

This shows that X„ i is a martingale difference array. Let’s now verify the conditions 
in McLeish’s theorem. 

1 Tl H” 1 

The variance is -- when a = b and 7 - 77 -— when a ^ b (see e.g. 

n + 2 {n-l){n + 2) ^ ^ 

n 

[ 2 ^ chapt. 4]). Thus, lim cr^ ^ = 1 and the array is already normalized. 

n—¥co ’ 


2=1 


By the row invariance of Haar measure, all the random variables X^^i are iden¬ 
tically distributed. Thus, Lyapunov’s condition for 6 = 2 simply states that 
lim rMyX^ J = 0. By Lemma EUl J = n^E[M^ = 0(1) and Lya- 


n—>-oo 

punov’s condition is satisfied. 


The last condition to show is that lim sup E E[X^ < 1. We split this into 




two sums: = E E Note that 

the random variables X^ iX^ j such that \i — j\> 1 are all identically distributed, 
and similarly for \i — j|= 1. The limit of the first sum is 


lim y 

n—yoo ' ^ 


E[Xl,Xl^] = lim E[rf^Xl^Xl^\ 

’ r/r, —’ ’ 


lim E[n'‘Mf 3 Mfc 24 ] = 1 . 

n—^oo 15 15 
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The asymptotic of the second sum is 

nxl,Xl^]=0{lln) 




and therefore 


limsup^E[X2,x2^^.] = 1. 




□ 


We’ve shown that if Cn is the matrix corresponding to the n-cycle (1, 2, 3,n) 
then each matrix entry {XICnXl'^)ab scaled by y/n converges to the standard normal 
distribution. It turns out that Cn is in fact a fairly generic permutation in ©„. The 
following lemma makes this precise. 

Lemma 2.3. Let M be an n x n Haar distributed orthogonal matrix and Z he a 
A/'(0,1) distributed random variable. Then 

lim ma,x\F{y/n{MQnM"’")ab < a;} — F{Z < a;}|= 0 

n—^oo Qji 

where the maximum is taken over the set of all n x n permutation matrices, Qn, 
such that the number of cycles in the corresponding permutation is at most 21ogn. 

Proof. Let d > 0 and let Qn be a permutation matrix whose corresponding per¬ 
mutation a G &n has kn < 21ogn cycles. For ease of notation, we set X„ = 
y/n{MQnM'^)ab and Yn = y/n{MCnM'^)ab- Since the distribution of Xn only 
depends on the cycle type of a, we can assume a = (1, 2, 3,..., ni)(ni -I- l,ni -I- 
2,..., n 2 )...(nfc„_i -I- + 2, ...,n) when decomposed into cycles. Then 

n n 2kn 

Xn — Yn = ^ ^ ^ Wi 

i—1 i—1 

where Wi = ^JnMa^aXIb^di for some indicies Ct and dt. For sufficiently large n, 
E[IFf] < Cjn for some constant C and by Cauchy-Schwarz, 

F[{Xn-Ynf]<A{kn)\C/n). 


By Chebyshev’s inequality, 

P{|X„ - y„|> e} < e-^F[{Xn - Ynf] < g-2 l6Clog^n 

n 

We also have the two inequalities 

^{Xn ^a;}<P{y„<a:-|-e}-|- P{|X„ — Fn,|> e} 

P{y„ < X - e} < F{Xn < x} +P{|X„ - r„|> e} 

Putting them together, we have 

nXn < x-e}-F{\Xn-Yn\> s] < P{X„ < x} < P{r„ < x + e} + F{\Xn-Yn\> e} 
Thus, 

|P{X„ < x} - P{r„ < x}| < P{F„ < X + e} - P{r„ <x-e} + 2F{\Xn - Yn\> e} 

< F{Yn < X -I- e} — P{Fti < X — e:} -|- C og n 


n 
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By the triangle inequality, 

^{Yn < a; + e} — P{F„ < a; — e} < < a: + e} — ¥{Z < a: + e}| 

+ ¥{Z <x + £} — ¥{Z < a; — e} + \F{Yn < x — e} — F{Z < x — £r}| 

Also, 

|P{X„ < x} - ¥{Z < x}|< |P{X„ < x} - P{r„ < x}| + |P{F„ < x} - ¥{Z < x}\. 

Since converges weakly to Z, |P{F„ < x} — P{Z < x}|< 5/2 for sufficiently 
large n. Choose small enough e such that P{Z < x + e} — P{Z < x — e} < 
i5/10. For large enough n, the three quantities |P{F„ < x + e} — P{Z < x + e}|, 

\¥{Yn < X — e} — F{Z < X — e}|, and ^ smaller than 5/10. 

n 

Then |P{Ai„ < x} — P{F„ < x}|< 5/2. Since Qn was arbitrary, this proves the 
lemma. □ 


Let Kn{a) denote the number of cycles in a random permutation a € ©„. The 


classical Goncharov theorem (see e.g. [1^ p. 116]) states that 

This implies in particular that P{iF„ > 21ogn} —?► 0. 

Now we are ready to prove: 


Kn - log n d 
Viogn 


A5(0,l). 


Theorem 2.2. Let M be a Haar distributed orthogonal n x n matrix and P a 
random nxn permutation matrix independent of M. Then 

Proof. Averaging over all permutation matrices, 

F{V^{MPM^)ab < x} = ^ y P{v^(MP,M^),fc < x} 

n\ 

For any 5 > 0 and sufficiently large n, Lemma 12.31 gives 
\F{^{MPM^)ab < x} - P{Z < x}| 

< A y \F{V^{MPaM^)ab < x} - P{Z < x}| 
n\ 

<5/2 + P{A:„ > 21ogn} 

<5 


□ 


3. Limiting joint distributions 

We can generalize Theorem l2.2l to obtain limiting joint distributions of the entries 
of MPM'^. By the Cramer-Wold Device, this is equivalent to determining the 
limiting distribution of arbitrary linear combinations of the matrix entries. In 
other words, for a sequence of n x n real matrices An, we wish to determine the 
asymptotic behavior of Tr{AnMPM"^). Theorem l2.2l savs that this trace converges 
to a standard normal when the (a, b) entry of A„ is ^/n and all the other entries 
are 0. 

We will now consider more general coefficient matrices A„. For convenience, the 
subscript will often be dropped. Normalize A so that Tr{AA’^) = n. 
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Remark 3.1. Recall that the limiting normality of Tr{AM) was proven in Q. A 
simple but key observation in the proof is that we can reduce to the case of di¬ 
agonal A by using the singular value decomposition. Writing A = UDV"^, we 
get Tr{AM) = Tr{UDV"^M) = Tr[DV^MU). This is equal in distribution to 
Tr{DM) by left and right Haar invariance. 

Unfortunately, we cannot reduce to the case of diagonal matrices in our situation 
since 

Tr{AMPM^) = Tr{UDV'^MPM^) = Tr{DV^MPM^U) 

= Tr{D{V'^ M)P{U^ M)^) 

and in general U ^ V. If A is normal, e.g. symmetric or orthogonal, then U = V 
and it would be possible to reduce to diagonal A. However, we will not make this 
assumption. 

As before, let C„ denote the n x n permutation matrix corresponding to the 
n-cycle (l,2,...,n). Then 

n n n 

Tr{AMCnM’^) = 

i —1 j — 1 k —1 

n n 

where Xn,k = It is easy to see that Xn,k is a martingale 

i=i i=i 

difference array. Squaring Xn,k and expanding, we get 


^n,k — ^2 ^2 ^ij^j,k^i,k+l 
i=l j=l 

n n n n 

+EEE AijAimMj^kMm,kM)^j._^_^ + EEE A,jAijMl^M,^k+iMi ,fc+l 

j^l m^j j^l 

n n 

+EEEE Aij AiYn.Mj^kMi^k+l^m,k^l,k+l 

i=l j=l l^i 

By Remark 12.11 the mixed moments are 

= 1/n^ + Oil/n^) 

E[Mj^k^i,k+iMi^k+i] = 0 

a j = i,m = I or j = l,i = "rn 
otherwise 


E[AIj,kM,,k+lMm,kMl^k+l] = < A-l)n{n+2) 


Therefore the variance is 


lE[^n,fc] - ( 72 + ( TI ) ) E E “ 


1 


i=i i=i 


(n — l)n{n -1- 2 ) 


'^2^{AijAji + AiiAjj) 




(3.1) 

Recall we have the constraint It is clear that is maximized 

when Aij = Aji and minimized when Ay = —Aji. In any case, the magnitude of 
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the sum is at most n. By the method of Lagrange multipliers, AaAjj achieves a 

maximum value of — n when A is the identity matrix and a minimum value of — n 
when ^ An = 0. Without any further restrictions, the sum AaAjj can fluctuate 

wildly. For there to be any hope of a limiting distribution for Tr{AMCnM’^), we 
- 7 ^ = c for some constant 0 < c < 1. Then we have the following: 

Theorem 3.1. Let An he a sequence ofnxn matrices such that Tr{A’^A) = n. As 
usual M denotes an n x n Haar distributed orthogonal matrix and Cn the permuta- 

tion matrix (1,2, If lim = c, t/ien TrlTMCnM^) —>■ VlO, 1 — c). 

n—foo ^^ 


Proof. From equation (EH), we see that the limiting variance of Tr{AMCnM"^) is 
given by 

n 





lim nW^Xn i 

n—^oo ’ 


= 1 — C 


Next, we show that Lyapunov’s condition for the martingale difference array 
Xn,k holds with <5 = 2, i.e. lim nE[Xn k[ = 0. 

n—>-oo ’ 

4 

When expanded out, X^ ^ contains terms of the form n ■^ipjp^jp,kMip,k+l- 

p=l 

Recall the only terms with non-zero expectation are those with an even number 
of matrix entries My in each row and column. Working through the possibilities, 
we get (using the shorthand i^j^lf^m to mean i,j,l,m all distinct and ^ to 
denote inequality up to an absolute constant as n —?► oo) 

E[<.]« E (,AaAjjAiiAnim. “t“ AijAjiAiiAnim AijAjiAiniAjnl 


■ Aij AjiAiiAnim AijAjiAiniAnii)^\AIi.kAIi,k+l-AIj^k-AIj^k+l-AIl^kAIl,k+lAlTn,k-AIm,k+l\ 


+ 'y ' + AaAjiAijAii -\- AijAjiAfi + AijAjiAnAu 

“b AaAjiAaAij 'jlK^AIi kAIi,k+i Mj^kMj^k+iMl^Ml^ 

,fc+l] 

+ y ^ (^AaAjjA^i + AaAijAiiAji + AjiAijAli)¥\Mi^kId[fj^j_iXIj^kMj^k+iMfj^ 

i=i^3 

jftm 

+ J2i^%AlnMlkMlk+iMl,k] 

3 
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By the discussion in Remark 12.11 it is easy to compute 

feM,-fe+i] = 0(1/71^) 

E[M,.fcM3fe+iM,-feM,.fe+iM^,] = 0{l/n^) 
E[M,.feM,,fc+iM3;,M3fc+J = 0{l/n^) 
E[M,.fcM3,+,iV4,iV4 fe+i] = Oil/n^) 

= 0{l/n^) 

= 0(l/n4) 

E[m4,m4,+i] = 0(l/n^) 

The sums involving the matrix entries of A are straightforward to bound. For 
example, we have 

AijAimAijAim < {A^jAij^ + A^^Aij)j2 < n 

j^m 

Going through all the terms, we see that E[X^= 0(l/n^) and Lyapunov’s 
condition is satisfied. 

The last condition of McLeish’s Theorem to show is 


limsup^E[X2_fcX2 ;] = (1 - cf. 

k^i 

(Note we didn’t normalize the array Xn,k)- By Cauchy-Schwarz, 

E[Xl^Xli] < E[x 4 _,] = 0 {l/n^). 

Thus, lim 1E[X^ f.X'^ J = 0. Therefore, we can assume \l — k\> 1. 

n—^oo • ^ ’ ’ 

ii-fei=i 

Consider the expectation of 




2 

n,fc 


X 


2 _ 

n,l 




■^ah^^h,k^a,k+l 


L a— 1 b —1 


/ n n Y 

EE AcdMdjMcj+\ j 

yc=i d=i J 


Again, going through all the possibilities, we have the following types of terms: 

( 1 ) ^^{AaaAdbAccAdd A AabAbaA^dAdc A AaaAbbAcdAdc)E\^XIa^kX[a,k+l^b,kX[b,k+l 

a^b 

c^d 

Mc,iAIc,i+iXId,iXId,i+i] 

( 2 ) 'y^XAcgAdbAcgAdb + AcaAdbAdgAcb)E[Mg^kXIc^k+iMb,kMd,k+iMajMc^i+iMbjMd,i+i] 

a^b 

c^d 

(3) E i^aa^cb^cb-^dd ~t“ ^aa-^cb^cd-^db “t“ Aca^ab^cd-^db)^[^a,k^a,k-\-l^b,k^c,k-\-l 

a,d^b,c 
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(■1) i: i^ca^db^ba^dc “t“ A.ca-^db^da-^bc H“ ^da-^cb^ba^dc 
a,d^b,c 

+ ^da^cbAda^bc)'^[Ma,kMc,k-\-lMty^k^d,k-\-l^a,l^b,l-\-lMc,lMd^l-\-l] 

( 5 ) 'Y2 i^aaAcbAdbAdc + AcaAabAdbAdc)^[Ma,kMa,k+lMb^kMc^k+lMbJMc^lMjl_^_J] 

d 

a^b^c 

( 6 ) ^ ] {AaaAbbA^^ + AabAbaA‘^^)¥,[Ma^kMa^k+lMb^kMb,k+lM^iMdi^i] 

c,d 

a^b 

(7) ^2i^caAcbAdaAdb)^[Ma^kMb^kM^i^_^iMa,lMb^lMdj+l] 

c,d 

a^b 

( 8 ) ^ {AlAlJE[Ml^Ml„+,MlMli^,] 

a,b,c,d 

The orthogonal integrals (using Remark 12.II) are given by 

E[Ma,kMa,k+iMb,kMb^k+iMcjMc^i+iMdjMdj+i] = 0(l/n®) 
E[Ma,kMc^k+iMb^kMd,k+iMa^iMc^i+iMb^iMd,i+i] = 0(l/n®) 
E[Ma,kMa,k+iMb^kMc^k+iMb^iMc^i+iMd^iMd^i+i\ = 0{l/n'^) 
E[Ma,kMc^k+iMb^kMd,k+iMa^iMb^i+iMc,iMd,i+i\ = 0{l/n'^) 

E[Ma,kMa,k+iMb^kMc,k+iMb^iMc,iMji_^_^] = 0(l/n®) 
E[Ma,kMa,k+iMb,kMb,k+iMliMli+^] = Oil/n^) 
E[Ma,kMb,kMl^^^Ma,iMb,iMli^^] = Oil/n^) 

E[Ml„Ml,,^,Ml,Mli+i] = 0{l/n^) 

Going through the 8 types of terms in the expansion of E[X'^ we see that 

the 2ncl, 3rd, 4th, 5th, and 7th terms are all 0(l/n^). The 1st is c^/n^ + 0{l/rA), 
the 6th is —cjn^ + 0(l/n^), and the 8th is 1/n^ + 0(l/n^). Note that in the 
expansion, we actually get two copies of terms of the 6th type since k and I can be 
switched. 

Thus, ElX^j^Xl i] = (1 - 2c + c^)/n^ + 0{l/n^). When |;-fc|> 1, Xlj^X^ are 
all identically distributed. This proves that limsup Y,E[Xl,Xl,] = {l-cf. □ 

k^i 

Theorem 13.II shows that Tr{AMCnM"'") A A/'(0, 1 — c) where Cn = (l,2,...,n). 
Unlike in Section[2l Cn is no longer a generic permutation. The asymptotic distribu¬ 
tion of Tr{AMQM"'") will depend on the number of fixed points of the permutation 
Q. The following lemma generalizes Lemma 12.31 from the previous section. 

Lemma 3.1. Let An be a sequence of n x n matriees such that Tr{A^A) = n and 



n—>-oo Ji 


for some constant s. Let M be an n x n Haar distributed orthogonal matrix and 

Z = j\r{o,i- s^). 
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Then for every non-negative integer f, 

lim max|P{Tr(AM(5„M^) < a:} — F{Z + /s < x}|= 0 

n—^oo 

where the maximum is taken over the set of all n x n permutation matrices, Qn, 
with f fixed points and at most 2 log n cycles. 

Proof. First, note that lim = s^. 

n^oo 

Let Qn be a permutation matrix with / fixed points and < 2 log n cy¬ 
cles. As before, we can assume the permutation has the form (1,2,..., ni)(ni -I- 
l,ni -I- 2,..., n 2 )...(nfc„_i -|- l,nfe^_i -I- 2,...,n). Set Xn = Tr{AMQnM'^) and 
Yn = Tr{AMCnM'^). Then the difference is 


where 


and 


2 kn—f f 

Xn-Yn= Y. I]C'4n.l+14n.2) 

m—1 m—1 


i=l j=l 

Vmp = AjjMj^kMj^k 


Vm,2 = Y^»^lk- 


Since E 




'■ i=l 


= 0(l/n), we have the convergence 


and therefore 


14 n ,2 ^ S 


Y ^ /«■ 


Let -I- Y^ Vm, 2 - For large n, we have the asymptotics 


nE[VY] =nYiK + A.,A,,)E[Ml^Ml^] = 0 ( 1 ) 


and 


nE[W^] = 0(1). 

Therefore, by Cauchy-Schwarz, 

E[iXn-Yff]<{2knf{C/n) 

for sufficiently large n for some constant C and therefore by Chebyshev, 
P{|X„ - Yf\> e} < e-\2knf{C/n) < 
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From the proof of Lemma 12.31 

< 4 - < a;}l < PK <x + e}- P{r^ < x - e} + 2P{|X„ - r^|> e} 

< F{Y^ < X + e} — F{Y^ <x — e} + ^ 

By the triangle inequality, 


P{y^ < a; + e} — F{Y^ < a; — e} < |P{F^ < x + e} — F{Z + fs < x + e}\ 

+ F{Z + fs < x + e} — F{Z + fs < x — e} + \F{Yf < x — e} — F{Z+ fs < a; — £}| 

By Theorem 13.11 (and the converging together lemma), Yf converges weakly to 
Z + fs. Thus, 

\F{Yf < x}-F{Z + fs < a:}|< S/2 

for (5 > 0 and sufficiently large n. Choose e small enough such that F{Z + fs < 
a; + £} — F{Z + fs < X — e} < (5/10. For large enough n, the three quantities 
\F{Y/ < a; + e} — F{Z + fs < x + e}|, \F{Yf < a; — e} — F{Z + fs < x — e}|, and 

^- 2 3 C og — bounded by (5/10. Then |P{X„ < a:} — F{Yl < a;}|< S/2 and 

n 

we have |P{X„ < x} — F{Z + fs < x}|< S for sufficiently large n. Since Qn was 
arbitrary, this proves the lemma. □ 


Putting everything together, we have the main result. 

Theorem 3.2. Let An, M, Z be as defined in Lemma \3.1l Let P be annxn random 
permutation matrix and let Y = Pois(/V) chosen independently from Z. Then in 
the limit of large n, 


Tr{AMPM^) A Z + sY. 


Proof. Let [/]„ be the subset of permutations in G„ with / fixed points and let 
Fn denote the number of fixed points of a random permutation in Recall that 

FuAy. 

Let (5 > 0. Choose a large integer F such that P(T > F) < 5/5. By Lemma IXTl 
for sufficiently large n, 


\F{Tr{AMPMA < x} - F{Z + sF < x}| 

^ oo oo 

iE E F{Tr{AMP^MA < x} - ^ P{F = f}F{Z + fs<x} 

/ = l(TG[/]„ /=i 


^ I 2 

< P{F„ >F}+ P{F > F} + Eb E F{Tr{AMP,jM^) < x} 

/=i ^ -e[/]n 



-p{r = /}P{F + /s<x} 
{F{Tr{AMP^MA < x} - F{Z + fs < x}) 

+ |P{F„ = /}-P{F = /}|^ 


(yC -I ^ I 

< — + — V V \F{Tr{AMP^M'^) < x} - F{Z + fs<x} 
3 n\ “ “ I 

/ = 1 (T€[f]n 
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The following corollaries illustrate a few special cases of this theorem. 

Corollary 3.1. Letk he a fixed positive integer. Let M be an nxn Haar-distributed 
orthogonal matrix and P be a random nxn permutation matrix. Then the joint 
distribution of k entries of the random matrix MPM"'" normalized by y/n is asymp¬ 
totically jointly i.i.d standard normal as n ^ oo. 

Corollary 3.2. Let An be a sequence of diagonal nxn matrices sueh that Tr{AX'") = 
n and An = y/lja for I < i < an for some parameter 0 < a < 1. Let 
Z = A/'(0,1 — a) and Y = Pois{l) be independent random variables. Then 

Tr{AMPM^) A Z + yf6tY. 

Remark 3.2. The result in Theorem l3.2l is for the defining representation of the sym¬ 
metric group. It would be of interest to see if similar distributional results for linear 
combinations of the matrix entries also hold for higher dimensional representations 
of the symmetric group. 


Acknowledgements: The author wishes to thank his PhD advisor Steve Evans 
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